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Abstract 

In order to solve an initial value problem by the variational iteration method, a 
sequence of functions is produced which converges to the solution under some suitable 
conditions. In the nonlinear case, after a few iterations the terms of the sequence 
become complicated, and therefore, computing a highly accurate solution would be 
difficult or even impossible. In this paper, for one-dimensional initial value problems, 
we propose a new approach which is based on approximating each term of the sequence 
by a piecewise linear function. Moreover, the convergence of the method is proved. 
Three illustrative examples are given to show the superiority of the proposed method 
over the classical variational iteration method. 
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1. Introduction 


The He’s variational iteration method (VIM) 0, 0 is a powerful mathematical 
technique to solve linear and nonlinear problems that can be implemented easily in 
practice. It has been successfully applied for solving various ODEs and PDEs 0, 0, 0 
iiQ. Convergence of the VIM has been investigated in some papers 0[i3. In [li], 
Zhao and Xiao applied the VIM for solving singular perturbation initial value problems 
and investigated its convergence. Recently, Wu and Baleanu in 0 have introduced a 
new way to dehne the Lagrange multipliers in the VIM for solving fractional differential 
equations with the Caputo derivatives. They also developed the VIM to g-fractional 
difference equations in [l^. The review article [l^ concerns some new applications of 
the VIM to numerical simulations of differential equations and fractional differential 
equations. 

Consider the following differential equation 


Cu{t) +N'u{t) = g{t), 


( 1 ) 
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where £ and M are, respectively, linear and nonlinear operators and g is an inhomo¬ 
geneous term. Given an initial guess UQ{t), the VIM to solve ([1]) takes the following 
form 

~\~ I t)(£'Uyi(s) -|- J\flimits') 5^(s))cis, Tfl 0, 1, 2, ... . (2) 

Jto 

where A is a general Lagrange multiplier (for more details see 0,08). 

As we see the VIM generates a sequence of functions that can converge to the 
solution of the problem under some suitable conditions. After a few iterations each 
term of this sequence often involves a dehnite integral whose integrand contains several 
nonlinear terms. In this case, computing a high accuracy solution would be difficult 
or even impossible by the common softwares such as MAPLE, MATHEMATICA or 
MATLAB. See 0 for an application of the VIM to quadratic Riccati differential equa¬ 
tions that illustrates this. To overcome this problem, Geng et ah in 0 have introduced 
a piecewise variational iteration method for the quadratic Riccati differential equation. 
Their main idea is that for solving the problem in a large interval, the interval is split 
into several subintervals and the problem is solved in each subinterval by the VIM 
in a progressive manner. In this way, the accuracy of the computed solution can be 
considerably improved. However, the main mentioned problem is still unsolved, since 
in each subinterval the VIM is directly implemented. 

In this paper, to solve one-dimensional initial value problems, we modify the VIM 
in such a way that each term of the sequence is interpolated by a piecewise linear 
function. Hereafter, the proposed method is called the IVIM (for interpolated VIM). 
In spite of the VIM, the IVIM does not need any symbolic computation and all of 
the computations are done numerically. Therefore, we can compute hundreds of the 
sequence terms in a small amount of time. 

The rest of the paper is organized as follows. The IVIM is introduced in Section |3l 
Gonvergence of the proposed method is discussed in Section [3l Section 0] is devoted to 
presenting three numerical examples. Some concluding remarks are given in Section 01 


2. The IVIM 


Gonsider the one-dimensional initial value problem 

f u'{t) = t e [a,T], 

\ u{a) = Ua- 


( 3 ) 


For simplicity, we assume that Ua = 0, otherwise one can use a simple change of variable 
it = u — Ua to obtain u{a) = 0. In this case, Eq. 02]) becomes 


Um+l{t)=Um{t)+ X{s,t) {u'^{s) - f{s,Um{s))) ds, m = 0, 1, 2, . . . , 


where Uo(t) satishes the initial condition Uo{a) = 0. By integration by parts, Eq. 
can be rewritten as ^ 

Gm(t^ / •) t)ds^ 


( 4 ) 

m 

( 5 ) 


2 


where 


G„i{t) = (1 + X(t,t))umit) - X{a,t)um{a), 


and 


dX 

TJ ('®) ^m('S) “1“ •^('5, t)/'(s,'Um(s)). 

OS 


In order to present the IVIM, we take a natural number n and discretize the interval 
[a, T] into n — 1 subintervals with step size h = {T — a)/{n — 1) and grid points 

ti = a + {i — l)h, i = 1, 2,..., n. 

Now, we dehne the well-known B-spline basis functions (see j^) of hrst-order on the 
nodal points t*, i.e. 


^i{t) = < 


^n{t) = 


t ti—\ 

ti—\ ^ ^ ti-, 

h ’ 

^2+1 ^ 
h ’ 

t 1 

0, 

t < L_1 OT t > t, 

t tn—1 

tn — 1 ^ t ^ tyi^ 

h ’ 

0, 

t ^ 15 


i = 2, 3,..., 77, — 1, 


on [a,T]. Each function tpi is equal to 1 at grid point ti, and equal to zero at other 
grid points. Let 

Xh = span{ 9 ?i : i = 2,3,... ,n}. 

Every in Xh is a piecewise linear function of the form 

n 

t e [a,T]. 

i=2 


Obviously, ctj is the value of at the grid point Xj, jm., aj = v’X^\ti) 


,(0 


(a) = 0, since (pi(a) = 0, i = 2,3,...,n. See |ll 


Furthermore, 
for more details about the 


properties of this kind of functions. 

We now return to Eq. ([5]). The aim is to recursively compute UiS. Let Um has been 
computed and we intend to compute Um+i- The main idea of our method is to compute 
a piecewise linear interpolation of Um+i in Xh, instead of computing Um+i{t). In 
fact, instead of Wm+i, we consider its projection onto Xh- To do so, we have 


= ^U^+i{tr)iPr{t), t e [a,T]. (6) 

r=2 


All we need here is to compute Um+i(ti), i = 2,3,... ,n. By using Eq. (j5]), we have 

^ t-ijds. (7) 

J a 


3 









Now, we can approximate the integral term in the latter equation by different methods 
of numerical integration. But, we again compute a piecewise linear interpolation of 
Hm{s,ti) in Xh as 


r=2 

Replacing by Hm\s,ti) in Eq. (j4]) yields 

n ti ri 


( 8 ) 


r=2 


r=2 


where firi = / (pr{t)dt, r = 2, 3,..., n. It is easy to see that 


"0, i < r — 1, 

h 

= < 2 ’ * = 

^ h, i > r + 1. 


Now, from this and Eq. ([8]) we get 


2—1 


^m +1 (^ 2 ) ~ (^ 2 ) Gffi (^ 2 ) h ^ Hyyi (tr-} ^ 2 ) ^ H^n (^25 ^ 2 )’ ^ 2,3,. (9) 


r=2 


Then, from ([6]) we set 

n 

ul^hit) ( 10 ) 

r=2 

Remark 1. Consider the special case that Cu = u' + au, for some « G M. This means 
that f{t,u(t)) = —au(t) —ATuit) + g(t), where J\fu(t) is the nonlinear term appeared 
in Eq. ([T]). In this case, we have \{s,t) = and therefore 1 + = 0. 

Obviously, from Eq. (|1]) we see that Um{a) = 0, for m = 0,1,.... Hence, Gmif) = 0, 
for m = 0,1,.... On the other hand, we have 

Hm{s,t) = e°‘^'‘~^\-aum{s) + f{s,Um{s))) 

= + aUmis) + Mum{s) - g{s)) 

= e^^^-^\Mu^{s)-g{s)). ( 11 ) 

3. Convergence of the method 

In this section, we present the convergence of the VIM and the IVIM. We hrst recall 
a lemma which gives conditions for the existence and uniqueness of a solution to the 
problem ([3]). 
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Lemma 1. Consider the initial-value problem where f : [a, T] x M —)■ M zs contin¬ 
uous with respect to the first variable and satisfies a Lipschitz condition with respect to 
the second variable, i.e., 


\f{t,u{t))-f{t,v{t))\<-i{t)\u{t)-v{t)\, te[a,T], (12) 


where ■y{t) is a nonnegative continuous function. Then, there exists a unique solution 
to the problem. 


Proof. See j^, 14 


To establish convergence of the VIM, we state the following theorem. 

Theorem 2. Let f : [a,T] x C^{[a,T]) ^ M be a continuous function with respect to 
its first variable which satisfies the Lipschitz condition TWt) . If Um(t) e C^([a, T]) for 
m = 0,1,2,..., then the sequence defined by converges to the solution of 

Proof By Eq. (|1]), we have 

Um+l{t) = f {{1 -I X{s,t))u'^{s) - X{s,t)f{s,U,n{s)))ds, 771 = 0,1,2,..., (13) 

J a 

since Um{o,) = 0. Let u(t) be the exact solution of Eq. (|3]). By Eq. ITT]) , it is readily 
seen that 

'^(t) = [ {{1 + X{s,t))u'{s) - X{s,t)f{s,u{s)))ds. (14) 

J a 

Let Emit) = Um(t) — u(t) be the error produced by the VIM at mth iteration. Since 
the hrst derivative of Eo(t) is continuous on the closed interval [a,T], there exists a 
constant M such that 


\E',{t)\ < M\Efit)\, Vte[a,T]. (15) 

From Eqs. flTdl) . flTTl) and flT^ . it follows that 

l^i(^)l = I [ {i^ +Hs,t))Eo{s) - X{s,t){f{s,Uo{s)) - f{s,u{s))))ds\ 

J a 

< j |1 + A(s,t)||Eo(s)|ds + f \'j{s)X{s,t)\\Eo{s)\ds. (16) 

J a J a 


Then, it turns out that by Eqs. flT5|) and ffT6|) . 


Efit)\<Ma 


|Eo(s)|cis, 


( 17 ) 


in which 


Ma = max{M max |l + A(s,t)|, max | 7 (s)A(s,t)|}. 

s,tG[a,T] s,t(^[a,T] 
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(18) 


Let llEjll = max \Ei(t)\, i = 0,1,.... Hence, 

te[a,T] 

\E,{t)\ < \\Eo\\M,{t-a). 

Therefore, similar to Eq. flTTI) and using Eq. ffTSj) . we deduce 

\E2{t)\ < Ma f \Ei{s)\ds < M^\\Eo\\ f {s - a)ds = \\Eo\ 
J a J a 

l^sWI <M, [ \E2{s)\ds < M^jEoW [ {s-afds = \\Eo 


' 2 ! 

' 3! 


and in general. 




ml " " ml 

Then \\Em\\ —)■ 0 as m tends to oo and the convergence of the method is established. 

□ 


Now, we are ready to prove the convergence of the IVIM. 


Theorem 3. Let f : [a,T] x C^{[a,T]) be two times continuously differentiable 
function with respect to the first variable and satisfies the Lipschitz condition fWfl . If 
the VIM of 0) is convergent, then so is the IVIM of /[lU) . 

n 

Proof. Dehne Umit) := J2umiti)ipi{t). For each grid point tr, we have Umifr) = 

i=2 

Um{tr). Also, by Eq. ([9]), we have 

f^m+l ) (1 T A (tj, tjffiu^YL (fi) h ^ ^ I s=tr (A) T ^ ifri tf) f (fr i ) 

r=2 

h d\ 

^ Qg l^=ii T f ifi^ ■ (^9) 

On the other hand, by Eq. (I5]) it is clear that the exact solution u{t) aX t = ti satishes 
the following relation 

u 

/ d\ 

(^(s, ti) u{s) + A(s, ti)f{s, u{s)))ds. (20) 

a 

Now, by the trapezoidal rule of integration, we have 



,d\ 


s,ti) u{s) + X{s,ti)f{s,u{s)))ds 


h dX 

2 <&' 

i—1 


S, ti) \s=ti u{ti) + X{ti, ti)f{ti, u{ti)) + A(a, ti)f{a, 0)) 


dX 


T/i ^ (s, ti) u{tfi) T A(fj., ti) f (t^, u(tr))) 


r=2 




12 ds2- 


( 21 ) 
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where ^ G {a,ti). Let EmiU) = UmiU) - u{ti) = Um{U) - u{ti) and EmiU) = Um{U) - 
u{ti). By Eqs. ([19]), (EOj) and (|2T]), we get 




2—1 


d\ 


1(1 + X{ti,ti))Ejn{ti) - h^{-^{s,ti) Em{tr) 


r=2 


h,d\ 


{tri'Um.iX'i'')') f ij'r t ^(j'r')')')') 2 ^ cls Ejji(tij 

+X{ti, Um{ti)) - f{ti, u{ti)))) - ^A(a, u)f{a, 0) 


12 ds2- 

Now, by the Lipschitz condition, it follows that 


d\ 

E^+i{ti)\ < |l + A(ti,L)| \E^{ti)\ +h'^{\ — {s,ti) I \Emitr)\ 

r=2 

h d\ 

+ \'y{tr)X{tr,ti)\ \Em{tr)\) + \s=ti \ \Em{ti)\ 

2 os 

+ b{ti)X{ti,ti)\ \Em{ti)\ + ||A(a,ti)/(a,0)| + “(^))l^=« 

^ El\Em{ti) \ + hK2 \Em{tr) \ + —K 2 \Em{ti)\ 

r=2 

+—\X{a,ti)f{a,0)\ H- -—h^K^, 


where |1 + A(t,t)| < iLi, ||j(s,t)| + | 7 (t)A(s,t)| < K 2 and \-^f{s,u{s))\ < K 3 for 
some constants Ki,K 2 and K^. Now, if the VIM is convergent, then \Em{ti)\ —)■ 0 as 
m ^ 00 and therefore, as a result when m —)■ cx) and h —)■ 0, we have |-E'm(L)| —t 0 for 
each grid point fj. It is necessary to mention that 


2—1 

lim lim h \ |Em(^r)| = hm 

m—^oo /i—>-0 ^ ^ m—>-cx) 

r=2 

Therefore the proof is completed. □ 

4. Illustrative examples 

In this section, we present three numerical examples to show the effectiveness of 
our method. In order to do so, we compare the results of the IVIM to those of the 
VIM. 




\Em{t)\dt = 0. 


'xi 


Example 1. We consider the quadratic Riccati differential equation of the form 


f u'(t) = 2 u(t) — v?{t) + 1, f e [0,1], 
\ m(0) = 0. 
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Table 1: CPU times (in seconds) for computing the approximate solutions by the VIM and the IVIM 
for Example 1. 


m 

1 

2 

3 

4 

5 

6 

7 

8 

VIM 

0.015 

0.093 

0.109 

0.204 

0.672 

2.485 

20.034 

383.469 

IVIM 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


The exact solution of fl2^ is given by (see for example [H) 

u(t) = 1 + y /2 tanh f \/2t H— log -] . 

\ 2 ^ V 2 + 1 J 

We assume that Cu = u' — 2u, Afu = v? and g = 1. By setting a = —2, according to 
Remark 1 we get Gm = 0 and 

H^{s,t) = e2(*“")(M^(s) - 1), 

for (s,f) G [0,1] X [0,1]. Therefore, Eq. (|5]) takes the following form 

Um+i{t) = [ -ul^{s))ds, 

Jo 

which is equivalent to the iteration produced by the VIM. In Figure [H the approximate 
solutions computed by the VIM and the IVIM for different values of m together with the 
exact solution are depicted. In the IVIM, we have considered n = 41. The hgure shows 
that the solution computed by the IVIM is in good agreement with that of the VIM and 
both of them converge to the exact solution. Table [U shows a comparison between the 
CPU times (in seconds) for computing the approximate solutions computed by these 
two methods. As we see, the CPU times for computing the approximated solution 
obtained by the VIM increase drastically as m increases, whereas those of the IVIM 
are negligible. 

For more investigation we set m = 10 and n = 33, 65,129, 257. In Figure[2]the log^^Q 
of the absolute error of the computed solutions by the IVIM, denoted by logj^Q Err, are 
displayed. The CPU times for all of the four values of n are 0.000. Similarly, for m = 40 
and n = 1000, 2000, 3000,4000, the results are shown in Figure |3l Here the CPU times 
for computing the approximate solutions by the IVIM are 0.265, 0.859,1.891, 3.203, 
respectively, for n = 1000, 2000, 3000,4000. These hgures show that the computed 
solutions obtained by the IVIM are in good agreement with the exact solution. More¬ 
over, the results show that the CPU times for setting up the IVIM is very small in 
comparison to the VIM. 

Example 2. We consider the initial value problem 


= Icost, t G [0,3], 
n(0) = 0. 


( 23 ) 








Figure 1: A comparison between the exact and the approximate solutions computed by the VIM and 
the IVIM for the Riccati differential equation for Example 1. 
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Figure 2: logj^p of the absolute error of the solutions provided by the IVIM for to = 10 and n 
33, 65,129, 257 in Example 1. 


m=40 



Figure 3: log^Q of the absolute error of the solutions provided by the IVlM for to = 10 and n 
1000, 2000, 3000,4000 in Example 1. 
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m=1 



Figure 4: Exact solution together with the approximate solutions computed by one iteration of the 
VIM and the I VIM for Example 2. 


The exact solution of the problem is u{t) = tv^sint. Let Cu = u', = 

— I cost and g = 0. Setting a = 0, from Remark 1 we get Gm = 0 and 

Hm{s,t) = COSS, 

for (s,f) G [0,3] X [0,3]. Therefore, from Eq. ([5]) the VIM takes the following form 


ttm+l (^) 



COSS ds. 


We have provided a MAPLE code for computing Wm’s- In Figure 0] we have depicted 
the exact solution together with the computed solutions by one iteration of the VIM 
and IVIM. As is seen, there is a good agreement between the computed solutions by 
the VIM and IVIM. For m > 2 the VIM could not provide an explicit expression for 
the Urn's, while the IVIM is well suited for sufficiently large values of m and n (Figures 
5 and 6). 

To see the efficiency of the IVIM, we consider the solutions provided by the IVIM 
with m = 10 and different values of n (= 33, 65,129 and 257). Figure [5] shows the log^g 
of the absolute error for the given solutions by IVIM. As we can see, the CPU times 
for computing these solutions are 0.000, 0.016, 0.031 and 0.047, respectively. Similar 
results have been shown in Figure E] for m = 40 and n = 1000,2000,3000 and 4000. 
The CPU times for computing the approximate solutions are 0.281,0.891,2.234 and 
3.484, respectively. 


Example 3. Although we have presented the IVIM for the hrst-order initial value 
problems, by this example we illustrate that the method can be implemented for the 
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Figure 5: log^g of the absolute error of the computed solutions by the IVIM for m = 10 and n 
33, 65,129, 257 for Example 2. 


m=40 



Figure 6: logj^g of the absolute error of the computed solutions by the IVIM for to = 40 and n 
1000, 2000, 3000,4000 for Example 2. 
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higher-order initial value problems. To do so, consider the second-order initial value 
problem 

( u"{t)-2u'{tf+ u'{t)+u{t) = g{t), te [0,1.5], . . 

\ m(0) = m'(0) = 0, ^ ’ 

where g{t) = t-|-2sin^ Ssin'^ |. The exact solution of this problem is u{t) = x — sint. 
Letting v{t) = u'{t) for t G [0,1.5], Eq. fl2T|) can be written as a system of two 
simultaneous hrst-order differential equations 


r u\t) = v{t), 

\ v'{t) - 2v{tf + v{t) -F u{t) = g{t), 


with the initial condition m( 0) = n(0) = 0. Having in mind the exact solution of 
(12TD . we see that the exact solution of the problem (1251) is given by {u(t),v(t)) = 
{t — sint, 1 — cost). 

To implement the VIM for solving fl25l) . in the hrst equation of fl25|) the linear and 
nonlinear terms are chosen as u' and —n, respectively. In the same way, for the second 
equation in fl25|) . v' + v and —2v‘^ + u are considered as the linear and nonlinear terms, 
respectively. Similar to the previous examples the VIM iteration takes the following 
form (for more details the reader can refer to [lil]]): 





(26) 


where uq and vq are two given functions satisfying mo( 0) = no(0) = 0. Obviously, from 
(126|) it follows that Um{0) = Vm{0) = 0. Therefore, by straightforward integration by 
parts Eq. (126|) can be written in the simpler form 


Um+l{t) = / V,n{s)ds, 

Jo 

Vm+l{t) = [ {‘2Vm{s)‘^ - Mm(s) + fi'(s)) ds 

Jo 

= [ e"~^2vm{s)‘^ + g{s)) ds - [ e"~^Umis)ds. 


To implement the IVIM, we set 


(27) 


Hm{s, t) = {2Vmis)‘^ + g{s)) , 

Hm{s,t) = 


In this case, Eq. (l27jl can be rewritten as 


Um+l(t) / Hm(^S,t)dSy 

Jo 

^m+i(t) / Hjji(^s,t)ds / , t)ds. 

Jo Jo 


(28) 
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Table 2: CPU times (in seconds) for computing the approximate solutions by the VIM and the IVIM 
for Example 3. 


m 

1 

2 

3 

4 

5 

VIM 

0.046 

0.313 

1.531 

16.750 

576.4 (Fail) 

IVIM 

0.010 

0.012 

0.012 

0.014 

0.015 


Taking Eq. (|5]) into account, one can approximate Um+i and Vm+i by the idea described 
in Section [3l Since the implementation of the method is straightforward, the details 
of the method are omitted here. Since, Hm involves the nonlinear term 2 n^, the main 
difficulties mentioned in Example 1 arise here, too. 

For the numerical results we use uo{t) = vo(t) = 0 for t e [0,1.5] as the starting 
point. We now present the numerical results of using the IVIM for solving the problem. 
To this end, we hrst take n = 20. In Figure [7] the graphs of Um and Vm for m = 
1,2, 3, 4, 5, computed by the VIM and the IVIM, together with u{t) = t — sint and 
v(t) = 1 — cost are depicted. This hgure illustrates that, in the IVIM, Um(t) — u(t) 
and Vm(t) —)■ v(t) for t E [0,1.5] as m tends to inhnity. Moreover, it shows that there 
is a good agreement between the solution provided by the IVIM and the VIM when 
the VIM gives a suitable solution. In addition, the VIM can not provides any suitable 
expression for ^ 4 , Um and Vm for m > 5. This is due to the nonlinear term involving 
Vm+i- For more investigation we report the CPU times (in seconds) for computing the 
solutions by the VIM and the IVIM methods in Table [2l As we can observe, the CPU 
times for the VIM are very small, whereas those of the VIM drastically grows as m 
increases. 

Finally, we set n = 1000 and compute the approximate solutions of the problem 
by the IVIM for m = 5 and m = 15. In Figure [H log^g of the absolute error of 
the computed solutions is displayed. It is noted that the VIM gives the approximate 
solutions in 0.082 and 0.201 seconds for m = 5 and m = 15, respectively. As we see, 
the IVIM provides quite suitable solutions for the problem in a reasonable amount of 
time. 

5. Conclusion 

We have proposed the interpolated VIM (IVIM) for solving one-dimensional initial 
value problems. The convergence of the method together with some numerical examples 
have been investigated. The numerical results show that the IVIM is more effective 
than the VIM. We have also shown the method can be easily implemented to the 
systems of ODEs. The following aspects can be considered in future works: 

1. Generalization of the method to two- or three-dimensional initial value problems; 

2. Implementation of the idea for the other methods such as the homotopy pertur¬ 
bation method; 

3. Using another type of basis functions such as the cubic spline functions and the 
Chebyshev polynomials. 
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m =1 m =1 



m=2 m=2 










Figure 7: A comparison between the graphs of u(t) = t — sint and v{t) = u'(t) 
respectively, with Um(t) and Vm(t), for m = 1, 2,3,4, 5. 


1 — cos(t) , 
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n=1000 


n=1000 




Figure 8: logj^Q of the absolute error of the solution computed by the IVIM for Example 3 for n = 1000 
and m = 5, 15. 


Acknowledgments 

The authors would like to thank the anonymous referees and the editors of the 
journal for their valuable comments and suggestions which substantially improved the 
quality of the paper. The work of the hrst author is partially supported by University 
of Guilan. 

References 

[1] S. Abbasbandy, A new application of He’s variational iteration method for 
quadratic Riccati differential equation by using Adomian’s polynomials, Journal 
of Computational and Applied Mathematics 207 (2007) 59-63. 

[2] D. Kincaid and W. Cheney, Numerical analysis, Brooks/Cole Publishing Company 
1991. 

[3] J.C. Butcher, Numerical methods for ordinary differential equations, 2nd edition, 
John Wiley and Sons, 2008. 

[4] F. Ceng, Y. Lin and M. Cui, A piecewise variational iteration method for Riccati 
differential equations, Comput. Math. Appl. 58 (2009) 2518-2522. 

[5] J.H. He, A new approach to nonlinear partial differential equations, Commun. 
Nonlinear Sci. Numer. Simul. 2 (1997) 230-235. 

[6] J.H. He, Some asymptotic methods for strongly nonlinear equations. International 
Journal of Modern Physics B 20 (2006) 1141-1199. 

[7] J.H. He, Variational iteration method - a kind of non-linear analytical technique: 
Some examples. Internal. J. Nonlinear Mech. 34 (1999) 699-708. 

[8] J.H. He, Variational iteration method - Some recent results and new interpreta¬ 
tions, J. Comput. Appl. Math. 207 (2007) 3-17. 


16 


















[9] J.H. He, G.C. Wu and F. Austin, Variational iteration method which should be 
followed, Nonlinear Science Letters A 1 (2010) 1-30. 

[10] J.H. He and X.H. Wu, Variational iteration method: New development and appli¬ 
cations, Comput. Math. Appl. 54 (2007) 881-894. 

[11] B. D. Reddy, Functional analysis and boundary-value problems: and introductory 
treatment, Longman Scientific & Technical, Harlow, 1986. 

[12] D.K. Salkuyeh, Convergence of the variational iteration method for solving linear 
systems of ODEs with constant coefficients, Comput. Math. Appl. 56 (2008) 2027- 
2033. 

[13] M. Tatar! and M. Dehghan, On the convergence of Hes variational iteration 
method, J. Comput. Appl. Math. 207 (2007) 121128. 

[14] A. Tavakoli and M. Ferronato, On existence-unigueness of the solution in a nonlin¬ 
ear Biot’s model. Applied mathematics and information sciences 7 (2013) 333-341. 

[15] C.C. Wu, New trends in the variational iteration method, Commun. Frac. Calc. 2 
(2011) 59-75. 

[16] G.C. Wu, and D. Baleanu, New applications of the variational iteration method - 
from differential eguations to q-fractional difference equations. Advances in Differ¬ 
ence Equations 21 (2013) 1-16. 

[17] G.C. Wu, and D. Baleanu, Variational iteration method for the Burgers’ flow with 
fractional derivativesNew Lagrange multipliers. Applied Mathematical Modelling 
37 (2013) 6183-6190. 

[18] Y. Zhao and A. Xiao, Variational iteration method for singular perturbation initial 
value problems. Computer Physics Communications 181 (2010) 947-956. 


17 


